Characterization of elastic scattering near a Feshbach resonance in rubidium 87 
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The s-wave scattering length for elastic collisions between 87 Rb atoms in the state \f, mj) = |1, 1) 
is measured in the vicinity of a Feshbach resonance near 1007 G. Experimentally, the scattering 
length is determined from the mean-field driven expansion of a Bose-Einstein condensate in a ho- 
mogeneous magnetic field. The scattering length is measured as a function of the magnetic field and 
agrees with the theoretical expectation. The width of the resonance is determined to be 0.20(3) G, 
the position of the zero crossing of the scattering length is found at 1007.60(3) G. 
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PACS numbers: 34.50.-s, 03.75.Nt, 32.80.Pj 



A Feshbach resonance offers a possibility for tuning 
the interactions in ultracold atomic gases simply by ap- 
plying a magnetic field. Feshbach resonances have been 
used to induce a controlled collapse of a Bose-Einstein 
condensate (BEC) 0], to create a coherent superposition 
of an atomic BEC and a molecular state 0, Lj , to realize 
a bright soliton in a BEC 0, 0], and to create a BEC 
of Cs atoms 0. Various groups are currently trying to 
use a Feshbach resonance to create a superfiuid phase in 
a degenerate Fermi gas. First results towards this goal 
have recently been published Q. 

Feshbach resonances have been observed in gases of 
various alkali atoms 0, IS S, IS E3, E3, E3, E3, E3, - 
With the isotope studied in this paper, 87 Rb, two tech- 
nical problems must be tackled. One problem is that in 
87 Rb Feshbach resonances exist only for internal states 
that cannot be held in magnetic traps (except for some 
spin mixtures). Another, more severe problem is that the 
Feshbach resonances in 87 Rb are unusually narrow: The 
broadest resonance located near 1007 G is predicted to be 
only 0.17 G wide. This imposes severe constraints on the 
current sources used to create the magnetic field. Yet, 
87 Rb is the isotope used in the vast majority of today's 
BEC experiments, thus stimulating a strong interest in 
investigating the possibility of manipulating the atomic 
interaction in this isotope. In a recent experiment, we 
were able to precisely locate more than 40 Feshbach res- 
onances in 87 Rb by observing enhanced atom loss |16| . 
Similar results for the broadest resonance were obtained 
very recently at the University of Oxford E3 ■ 

In the ultracold regime, the elastic scattering prop- 
erties are fully characterized by the s-wave scattering 
length a. Quantitative measurements of a near Feshbach 
resonances have previously been reported for 3 isotopes: 
23 Na HE3, 85 Rb USES 13, and 40 K [l^. These exper- 
iments are based on a measurement of either the mean- 
field energy of a BEC or the thermalization rate in a 
non-degenerate gas. 

In this paper, the mean-field energy of a BEC is used 
to measure a in 87 Rb in the vicinity of the Feshbach res- 
onance near 1007 G. Our method differs slightly from 
the one previously published |SL llSL Il9| . There, the final 
value of the magnetic field was applied while the atoms 
were still trapped and the atoms had time to equilibrate 



before being released from the trap. In this paper, how- 
ever, we first turn off the trap and then switch to the final 
value of the magnetic field. This improves the signal-to- 
noise ratio as discussed at the end of this paper. 

The experimental set-up is only briefly summarized 
here. More details are given in Ref. LJ. A double- 
MOT system is used for cooling and trapping of atomic 
87 Rb. The atoms are then optically pumped to the state 
|/, rrif) = |1, — 1) and loaded into a Ioffe-Pritchard mag- 
netic trap, in which a radio-frequency (rf) field drives 
evaporative cooling to the BEC phase transition. It 
turned out that BECs with the same properties as in 
Ref. [3 

can be created with a much faster rf sweep than 
described there. The optimized rf sweep lasts only 5.1 s, 
resulting in an overall cycle time of 15 s. 

After creation of the BEC, the atoms are transferred 
into an optical dipole trap and held there for roughly 
0.5 s. The optical trap is made of two beams from a 
Nd:YAG laser (A = 1064 nm), with the beams crossing 
at right angles. One beam propagates horizontally, along 
the symmetry axis of the magnetic trap. The second 
beam subtends an angle of 25° with the horizontal plane. 
The beam waist (l/e 2 -radius of intensity) and power are 
33 /im and 38 mW for the horizontal beam; and 77 /im 
and 115 mW for the second beam. Due to the tighter 
waist, the horizontal beam creates a stronger confinement 
than the second beam. 

The horizontal beam creates an estimated trap depth 
in the horizontal plane of ~ ks X 3 /iK. In the verti- 
cal direction, the trap barely supports the atoms against 
gravity with a trap depth of ~ ks x 0.8 /iK. The horizon- 
tal beam creates negligible confinement along its propa- 
gation direction. The second beam solves this problem 
without changing the confinement in the other two di- 
rections significantly. The polarizations of the two light 
beams are orthogonal, so that no optical lattice is formed. 

In a coordinate system, where the z axis is the sym- 
metry axis of the magnetic trap and where gravity points 
along x, the trap frequencies are (u> x , uj v ,uj z ) — 2ttx (120, 
170, 50) Hz. The first two frequencies differ due to the 
gravitational sag. Typically, the atom numbers in the 
BEC and in the thermal fraction are 10 5 each. The peak 
density in the BEC is - 2 x 10 14 cm~ 3 . The optical 
trap is not designed to maximize the atom number in 
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the BEC. Instead, we try to avoid high peak densities in 
order to reduce density-dependent losses. 

During the transfer from the magnetic trap to the opti- 
cal trap, a small magnetic bias field (B « 1 G) preserves 
the spin polarization of the atoms. After the transfer 
to the dipole trap, a strong homogeneous magnetic field 
(B w 1000 G) is applied in the opposite direction. The 
rapid switching-on of this field (at rates up to 700 G/ms) 
flips the atomic spins from state 1) to the absolute 
ground state |1, 1), in which the Feshbach resonance oc- 
curs. With this method, the spin-flip efficiency is so close 
to 100 %, that a Stern-Gerlach method does not show any 
significant signal from atoms in wrong spin states. We 
estimate the detection limit of this measurement to be 

The 1000-G field is typically not set exactly onto the 
Feshbach resonance right away. Instead, the field is held 
a few G above or below the Feshbach resonance for typi- 
cally 500 ms. During this time, thermal drifts due to the 
^12 kW heat load dissipated in the coils have some time 
to settle. Next, the optical trap is switched off and simul- 
taneously B is jumped to a value very close to or right 
at the Feshbach resonance. After holding B at its final 
value for 6 ms, B is switched off completely. The jump to 
the final field value actually takes ~ 0.5 ms, the complete 
turn-off ~ 2 ms. After an expansion time of the order 
of 20 ms, an absorption image of the expanded cloud is 
taken with a CCD camera. The atom number and size 
of the BEC and of the surrounding thermal cloud are de- 
termined from a two-dimensional fit to the CCD picture. 
The magnetic field was calibrated using microwave spec- 
troscopy in the vicinity of the Feshbach resonance with 
an accuracy of 0.03 G. 

Figure ^ shows the size and atom number of the ex- 
panded BEC measured with the timing sequence de- 
scribed above. As expected, the size increases (decreases) 
as the Feshbach resonance is approached from below 
(above). Additionally, the atom number is reduced as 
one approaches the resonance from either side. No BEC 
is left between 1007.37 G and 1007.53 G. The separation 
of adjacent data points in Fig.QJis 13 mG which roughly 
equals the estimated magnetic field resolution due to the 
measured current noise. 

As one gets very close to the Feshbach resonance from 
either side, the upward (downward) trend in the width is 
reversed. On the low-field side of the resonance, this is a 
trivial consequence of the decrease in atom number, be- 
cause a lower-density BEC releases less mean-field energy 
that drives the expansion. On the high-field side of the 
resonance, however, the effect is related to the instability 
of the BEC in a regime of negative a discussed later in 
this paper. 

In order to extract a from the measured size of the 
expanded BEC, the expansion process must be modelled. 
To this end, the model presented in Ref. [2l| is extended 
to include a possible time dependence of a during the 
expansion. We start by first ignoring the observed change 
in atom number. This effect will be incorporated in the 
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FIG. 1: Vertical half width (•) and atom number (o) of the 
expanded BEC as a function of magnetic field B. Data points 
in the left (right) half of the figure were obtained by jumping 
towards the Feshbach resonance from smaller (larger) B, and 
for an expansion time of 18 ms (23 ms). B was jumped to 
its final value at the moment of release of the atoms from the 
trap. B was held there for 6 ms and then switched off. The 
vertical line at 1007.60 G indicates the onset of instability of 
the BEC due to negative values of a. 



analysis later on. 

Initially, the BEC is confined in a harmonic potential 
with trap frequencies (u> x ,u) y ,u) z ). In the Thomas-Fermi 
approximation to the Gross-Pitaevskii equation, the ini- 
tial density distribution n(r) is an inverted parabola 



n(r) = no 



E 

k=l 



(1) 



and n(r) = if the above expression is negative. Here, 
n = 15 N/(8nW x W y W z ) is the peak density. The half 
widths along the coordinate axes are 
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where m is the atomic mass, N the atom number in the 
BEC, and a* the initial scattering length before release. 
At time t = 0, the trap is switched off, and a possible 
time dependence a(t) may begin. 

According to Ref. |2l|, the BEC profile stays a para- 
bola during the expansion. Its widths are scaled by the 
parameters 



Afc(t) 



W k (t) 
W k {0) 



(3) 



The evolution of the scaling parameters A^ is described 
by a set of coupled differential equations, Eq. (11) in 
Ref. [2l| • One can easily show that a time dependence of 
a during the expansion modifies this equation to 



a(t) 



1 
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The initial conditions are Afc(O) = 1 and Afc(O) = 0. a(t) 
appears in Eq. (J3j, because the expansion is driven by the 
mean- field energy which is proportional to a(t). cti ap- 
pears in Eq. because a% determines the initial widths 
Wk(0) with respect to which the A& are defined. 

As mentioned earlier, jumping a from its initial to its 
final value aj takes ~ 0.5 ms. We include this delay in 
the model in the following way 



a(t) = 



a f 

CL; - 



0.5 ms < t < 6 ms 
otherwise . 



(5) 



Here, ab g is the background value of the scattering length 
far away from the resonance and t = is defined by the 
turn-off of the optical trap. The fact that a is switched 
back to <2b g at t = 6 ms has little effect, because for 
the parameters of the experiment almost all mean-field 
energy is released during the first 6 ms (except for very 
small values of af). 

The atom loss visible in Fig. must also be included 
in the model. This is fairly easy if it is possible make the 
approximation that the atom loss is independent of the 
atomic density. With this approximation, the shape of 
the BEC remains parabolic during the expansion. With 
a time dependence of the atom number N(t), Eq. (jJJ is 
then modified to 



A fc = 



a{t)N(t) 



1 



Q>iNi A x Aj,A 2 Afc 



(0) 



Additionally, in Eq. |J2}, N is replaced by the initial value 
Ni . Since the dominant loss mechanism is actually likely 
to be density dependent, this approximation will only be 
reasonable if a small fraction of the atoms is lost. 

It is clear from Eq. JfjJ that the time dependence of the 
atom number is crucial. By varying the hold time at the 
final magnetic field from 1 to 6 ms, we experimentally 
checked for such a time dependence at B = 1007.35 G 
(last non-vanishing data point on the low-field side in 
Fig-D)' but found none: The final atom number was in- 
dependent of the hold time over this range. It is hence 
obvious that the loss does not occur continuously during 
the total hold time. Instead, the loss occurs during the 
first 1 ms of hold time. The origin of this loss might be re- 
lated to the formation of molecules, but this is not clear. 
Similar but weaker loss on short timescales was observed 
for a thermal cloud of atoms in our previous experiment 
16]. For simplicity, we assume that the atom loss occurs 
instantaneously at the same time as the change in a 
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0.5 ms < t 
otherwise . 
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For processing the data, the initial atom number Ni 
and the background value of the scattering length db g 
are needed. Ni was determined from the overall absorp- 
tion in a measurement far away from the Feshbach res- 
onance. In order to determine <Zb g , the size of an ex- 
panded BEC far away from the Feshbach resonance was 
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FIG. 2: Scattering length as a function of magnetic field. 
The solid line is a fit to the experimental data (•). The best 
fit yields AB = 0.20 G for the width and B pcak = 1007.40 G 
for the position of the resonance. Some experimental data 
points (o) were not included in the fit, because they lie in the 
regime of negative a where the BEC is unstable. 



measured. Using Eqs. with a(t) = a* = ab g , the 

background scattering length can be extracted, yielding 
abg = 108(30) a , where a is the Bohr radius. This is 
consistent with the theoretical value ab g = 100.5 ao for 
the state |1, 1} 22]. For further data processing, the the- 
oretical value was used. 

Given Ni and Ob g , W x (0) is calculated using Eq. J2J, 
and thus the measured expanded widths W x (t) can be 
converted into X x (t) according to Eq. 10 . For every data 
point, the observed final atom number Nf in the BEC 
was then used to numerically solve the coupled differen- 
tial equations JfjJ) for the experimentally applied sequence 
Eqs. 0. In the calculation, a/ was varied until the 
experimentally observed width was matched. 

With this method, the data shown in Fig^were pro- 
cessed to extract the scattering length, which is shown 
in Fig. [21 Figure [3 also displays a fit to the theoretical 
expectation (see e. g. [23)) 



a bg 1 



AB 



B-B 



peak 
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Here, AB is the width of the Feshbach resonance and 
Speak the position of the pole in a. 

Some data points (o in Fig.|2J) are not included in the fit 
for the following reason: Approaching the resonance from 
above in Fig.^ one can see that the size of the expanded 
BEC decreases as expected, until at 1007.60 G (vertical 
line in Fig. ^) this trend is reversed. At the same field, 
significant atom loss suddenly begins. We interpret this 
as the onset of instabilit y of the BEC in a regime of nega- 
tive scattering length [Il2 ll25ll26ll27l l28 | . This interpre- 
tation is further supported by the fact that the extracted 
value of a shown in Fig. [21 reaches zero at this field. For 
the parameters of our experiment, the critical scattering 
length for instability 0i|23 is a crit ~ — 10~ 3 ab g which is 
almost zero. Hence, we conclude that the zero crossing of 
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a is located at B zcro = B pcak + AB = 1007.60(3) G and 
we use only one free parameter AB in our fit. A more 
detailed analysis of the region of instability is clearly be- 
yond the scope of this paper. 

The fit in Fig. [21 is in reasonable agreement with the 
experimental data. Data points far away from the reso- 
nance are systematically too close to unity. This might 
be due to the fact that the magnetic field was not held 
infinitely far away from the resonance before jumping 
to the resonance. The best-fit value for the width is 
AB = 0.20(3) G, resulting in £ pcak = 1007.40(4) G. 
These values are consistent with the theoretical predic- 
tions AB theory = 0.17(3) G and B^°° ry = 1008.5(1.6) G 
[ill. Speak is also near the peak of the atom loss at 
1007.34(3) G measured in Ref. 0. 

As mentioned in the introduction, the method to de- 
termine a described in this paper differs slightly from the 
method published previously H, 0,0]- There, the final 
magnetic field was applied while the BEC was still in the 
trap and the system had time to equilibrate before the 
BEC was released from the trap. No matter if the expan- 
sion time is almost zero |l9| or long |3, [l8| (in which case 
B was still on during the initial expansion) , the observed 



BEC size W yields a oc W 5 . With the method used in 
this paper, however, a oc W 2 . This is because here the 
initial BEC peak density n, is independent of B. Hence, 
the mean-field energy (oc rija) is converted into kinetic 
energy, so that the final velocity (and thus W) is pro- 
portional to a 1 / 2 . Therefore, with our method noise in 
the determination of W is not amplified as much when 
extracting a. A similar technique was used in Ref. |(|, 
but there no quantitative values of a were extracted. 

To summarize, the mean-field driven expansion of a 
BEC was used to measure the scattering length a in 87 Rb 
as a function of the magnetic field in the vicinity of the 
Feshbach resonance near 1007 G. The position and width 
of the resonance were extracted from the data. Despite 
the fact that the resonance is very narrow, we have clearly 
demonstrated the variation of a over a wide range. The 
range explored here is currently not limited by magnetic 
field noise, but rather by fast atom loss near the reso- 
nance. 
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